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■ Abstract. The intricate interplay of atmospheric shock waves and a complex, 
variable radiation field with non-equilibrium dust formation presents a consid- 

, erable challenge to self-consistent modelling of atmospheres and winds of AGB 

|"T i ■ stars. Nevertheless it is clear that realistic models predicting mass loss rates and 

synthetic spectra are crucial for our understanding of this important phase of 
stellar evolution. While a number of questions are still open, significant progress 
has been achieved in recent years. In particular, self-consistent models for atmo- 
spheres and winds of C-stars have reached a level of sophistication which allows 
direct quantitative comparison with observations. In the case of stars with C/O 

■ < 1, however, recent work points to serious problems with the dust-driven wind 
■^j- | scenario. This contribution analyzes the basic ingredients of this scenario with 
■^j- , analytical estimates, focusing on dust formation, non-grey effects, and diffcr- 
PsJ ' ences between C-rich and O-rich environments, as well as discussing the status 
C of detailed dynamical wind models and current trends in this field. 

o 
6 

' Mass loss is one of the most pronounced features of AGB stars, influencing both 

their evolution and their observable properties in a decisive way. The most 
commonly accepted scenario is that of a pulsation-enhanced dust-driven wind: 
stellar pulsation causes atmospheric shock waves which intermittently lift gas 
above the stellar surface, creating dense, cool layers where solid particles form. 
The dust is accelerated away from the star by radiation pressure, dragging the 
gas along. The composition of the grains is determined by the relative elemental 
abundances in the atmosphere, with C-rich stars forming mostly amorphous 
carbon grains, and O-rich objects producing silicate dust. 

Current knowledge suggests that winds of C-rich AGB stars are rather 
well understood. For their numerous O-rich counterparts, however, recent work 
demonstrates that the opacities of silicate grains are insufficient for driving winds 
(Woitke 2006a and this volume), contrary to previous expectations. This prob- 
lem becomes apparent with the introduction of frequency-dependent radiative 
transfer in time-dependent models with a detailed description of dust formation. 
Ironically, this important step towards more realistic models which brought C- 
rich models into good agreement with observations, causes serious problems with 
the wind mechanism in O-rich stars. 

The core of this issue can be understood with comparatively simple analyti- 
cal arguments. A significant part of this review is devoted to isolating the crucial 
ingredients and analyzing them with simplified analytical descriptions, in par- 
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ticular non-equilibrium dust formation, with special attention to the influence of 
non-grey effects and differences between C-rich and O-rich stars. These simple 
estimates are compared to the results of existing numerical models, rounding 
up with an overview of the present status of such detailed models, and current 
trends. 

Detailed comparisons of the input physics of various numerical models and 
summaries of the historical developments in this field can be found in earlier 
reviews, e.g., by Willson (2000), Woitke (2003) and Hofner (2005). 

2. Basic ingredients and characteristic numbers 

Dust grains forming in the extended atmospheres of AGB stars condense under 
non-equilibrium conditions with temperature acting as a threshold, prevailing 
densities and abundances determining the efficiency of grain growth, and pul- 
sation and shock waves setting the time scales. In other words, condensation 
will usually be incomplete, and the dust-to-gas ratio is not a simple function of 
abundances. 

2.1. Grain temperature and condensation radius 

A prerequisite for condensation or survival of existing grains are temperatures 
below the stability limit of a specific condensate. In principle, a distinction 
between the respective temperatures of ambient gas and grains has to be made. 
In practice, however, gas temperatures tend to be lower than grain temperatures 
at a given point, and therefore the limit will usually be defined by the grain 
temperatures which, in turn, are given by radiative equilibrium, i.e. 



Using the assumption of a Planckian radiation field (with stellar surface tem- 
perature T* at J?*), geometrically diluted according to distance from the stellar 
surface 



in the equation of radiative equilibrium, the grain temperature as a function 
of distance from the star can be estimated if the wavelength-dependent opacity 
of the dust particles k\ is known. The latter can often be approximated in the 
critical wavelength region around the stellar flux maximum by K\ oc A~ p where 
the exponent p is characteristic of the grain material, leading to 








/ °° K\B x (T d )d\ 
f™K X B x (T*)d\ 




4 ..Planck 





Planck 



(3) 



where the superscript 'Planck' denotes a Planck mean of the opacity (see, e.g., 
Lamers & Cassinelli 1999 for details). The condensation radius r c , defined as the 
point where the grain temperature Td is equal to the condensation temperature 
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T c (stability limit) of the respective grain material, is therefore given by 



where we have used the approximation W(r) ~ (i?*/2r) 2 for r S> R*. 

Therefore, to determine the location where condensation may start for a 
particular material, we need to know its condensation temperature and the de- 
pendence of the opacity on wavelength around the flux maximum. For amor- 
phous carbon grains we have T c ~ 1500 K and p ~ 1. Assuming a stellar surface 
temperature of T* ~ 3000 K we obtain r c /R* ~ 3 which is in good agreement 
with detailed frequency-dependent models where dust usually forms at 2 — 3 i?* . 

For silicate grains the picture is more complicated. To illustrate the prob- 
lem, we consider olivine, i.e. a material with the composition Mg2 :r Fe2(i_. r )Si04 
(0 < x < 1) and a condensation temperature of T c ~ 1000 K in the relevant 
density range. On the one hand, the lower value of Tj/T* ~ 1/3 (compared 
to « 1/2 for carbon grains) leads to a stronger dependence of r c j R* on p (see 
Fig- 1])- On the other hand, the value of p varies strongly with the relative Fe- 
content (x) of the material, with larger (more positive) values of p for iron-rich 
grains. For forsterite (MgSiO^, representing one end of the chemical spectrum, 
lab data indicates p ~ —1 in the relevant wavelength rage (see Fig. 2 in Andersen, 
this volume). For MgFeSi04 (with equal amounts of Mg and Fe, corresponding 
roughly to solar composition) the value is closer to p ~ 2. Therefore forsterite 
may condense at a distance of ~ 3R* while MgFeSi04 will probably not form 
closer to the star than m 14 R*. Already these first estimates lead us to expect 
that silicate grains forming in AGB stars will tend to be iron-poor, a conclusion 
which is supported by other arguments and detailed modelling (see below) J3 

2.2. Dust opacity, shock waves and wind regimes 

Assuming for the moment that a certain dust species will actually form beyond 
its condensation radius, will the grains significantly contribute to driving a wind? 
This is, of course, a matter of the corresponding opacity produced by the grains. 
If we can assume that gravity and radiation pressure are the only relevant forces 
beyond r c , the (co- moving) equation of motion for a matter element can be 
written as 



where r is the distance from the stellar center, u the velocity, go = GM^/tq 
denotes the gravitational acceleration at the (arbitrary) point ro, kh the flux 
mean opacity, and M* and L* are the stellar mass and luminosity, respectively 
(all constants have their usual meaning). The material properties (and relative 
abundances) of the dust species determine the value of T which, in turn, de- 
termines the further dynamics of the matter element. Several special cases are 
immediately apparent: 



1 Similar arguments hold for pyroxenes, another possibly abundant grain material, with a com- 
position Mg^Fefi.^jSiOa (0 < x < 1). 
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Figure 1. Condensation radius (in units of stellar radius) as a function of 
p. The lower curve corresponds to T c /T* = 1/2 (representative of carbon 
grains; full line) and the upper curve to T c /T* = 1/3 (characteristic of sili- 
cates; dashed). The symbols show the approximate location of various mate- 
rials in this diagram (diamond: amorphous carbon; circle: Mg2Si04; square: 
MgFeSi04). Models with grey radiative transfer (or grey dust opacities) cor- 
respond to p = 0. 



• If no (or very little) dust is forming, T will be (close to) zero. The matter 
element will follow a ballistic trajectory, according to its initial velocity^ 
reaching a maximum distance r max of 



= 1 - ( ) where u esc = 

'"max \ U csc / 

and then falling back towards the stellar surface. 

• If T = 1, the r.h.s of the equation of motion vanishes. The matter element 
continues to move at a constant velocity. 

• If T > 1, du/dt is positive, the matter is accelerated away from the star. 

• If < T < 1, the final fate of the matter element (escape or fall-back) 
depends on the velocity at which it is moving when it reaches r c . 

The dynamical behavior in the regime < T < 1 can be investigated with 
the following simplified model: if we assume that T is constant outside r c (and 



^ (6) 



2 In the discussion here we always assume that the velocity before reaching the dust condensation 
radius is below the escape velocity. Otherwise it would be misleading to talk about a dust- 
driven wind. 
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zero inside), the factor (1 — T) is a constant that can be multiplied with go, 
resulting in an equation of motion 

Tt = ~ 9T ^f) 2 With 9Y = 9o{l-T) (7) 

which formally looks like the ballistic equation of motion for T = 0, but with 
a re-scaled gravitational acceleration gr, and therefore a re-scaled (lower) es- 
cape velocity u r sc = ^2GM*(\ — Y)/r. Consequently, the actual velocity of the 
matter element when reaching r c may be above u^ sc (while still being below the 
escape velocity without dust), resulting in a wind. 

To find the critical case that divides fall-back from outflow, we assume that 
the velocity prior to condensation is only due to the pulsation-induced shocks 
propagating through the atmosphere. We use the simplified picture that a shock 
passing through the matter element under consideration can be treated as an 
instantaneous acceleration to a velocity uq at the point ro (which will in practice 
be close to i?*), followed by a ballistic movement of the shocked gas. We consider 
the solution of the ballistic equation of motion (r = 0) for r < r c , requiring that 
the velocity at the point r c matches Ug SC , i.e., 

ul-u\r c ) = -2GM if (---) and u 2 (r c ) = (1 — T) . (8) 

For < T < 1 we can therefore distinguish two cases: 

• a pulsation-supported dust-driven wind regime for 

1>r >!i(l-(^) 2 ) or !>i* > Jl-rfi) (9) 

ro \ \UcscJ J Uesc V ^ rc ^ 

• a 'parachute regime' (decelerated fall-back) for 

-fl-f— ) 2 )>T>0 or J 1 - T ( r J>)>«L>0 (10) 

TO \ \UcscJ J V V r c/ U esc 

Assuming that ro ~ R* and, consequently, uo/u csc ~ ("U s hock/^esc) r= R the con- 
dition for a pulsation-supported dust-driven outflow can be reformulated in the 
more convenient form 




where u s h oc k denotes the velocity of the shocked gas (not to be confused with 
the shock amplitude). 

In order to determine the value of V for various grain materials, we first 
need to calculate the corresponding (flux mean) opacities. For simplicity, we 
assume that kh can be approximated with the small particle limit opacity of 
spherical grains at A max ~ 1/xm (flux maximum of the star). Therefore, we have 

K H ~ - / a 2 vrQ cx t(Amax) n(a)da = -Q ext (A ma x) / a 3 n(a) da (12) 
p Jo p Jo 
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where a denotes the grain radius, n(a) the number density of grains and Q ex t 
is the extinction efficiency. For grains small compared to the wavelength, the 
quantity Q' cxt = Q cx t/a becomes independent of the grain radius and can there- 
fore be taken out of the integral. The last integral represents the fraction of a 
given volume which is occupied by the grains, apart from a factor 47r/3, and 
it can be rewritten in terms of the space occupied by a monomer (basic build- 
ing block) in the condensed material times the number of monomers found in a 
certain volume, 

/ a 3 n(a) da = ^- V mon n mon = ^- Amonmp f c £c nR (13) 

J 47T 47T Pgrain 

where we have expressed the monomer volume Vmon in terms of the atomic 
weight of the monomer A mon and the density of the grain material p gr ain, and the 
number of monomers in a volume n mon by the abundance of the key element of 
the condensate e c , the degree of condensation of this key element f c , and the total 
number density of H atoms nn (m p = proton mass). Using bh = p/fl+Asue) mp, 
we finally obtain 

K H S3 0.5 Qcxt(^max) fc^c- (14) 

Pgrain 

Below, we list the properties of different grain materials and the corresponding V, 
assuming L/M = 5000 Lq/Mq and a degree of condensation of the key element 
f c = 1 (note that the value of T scales linearly with each of these two factors). 
The abundance of the key element e c corresponds to either all C not bound in 
CO for carbon grains (C/O = 1.5), or to Si for silicates. The data for Q' ext is 
taken from Rouleau & Martin (1991, amorphous carbon) and Jager et al. (2003, 
silicates; cf. Figs. 1 and 2 in Andersen, this volume). 

material -^mon Pgrain Qcxt £ c ^L/M=5000 

[g/cm 3 ] [1/cm] 
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amorphous carbon 12 1.85 2 • 10 4 3.3 • 10 4 
Mg 2 Si0 4 140 3.27 2 • 10 1 3.6 • lO" 5 6 • 10~ 3 

MgFeSiC-4 172 3.71 7 • 10 3 3.6 • 10~ 5 2 



As expected, the value of T for amorphous carbon grains is well above the 
threshold for a dust-driven wind, even for a degree of condensation f c ~ 0.3 — 0.5 
as typically found in detailed models. For silicate grains, the picture is, again, 
more complex: for forsterite (Mg2Si04), on the one hand, T is so far below the 
critical value that even the scenario of a pulsation-supported dust-driven outflow 
as discussed above is unrealistic, since it would require (« s hock/'"esc) r= R < , close 
to unity. Silicate grains containing equal amounts of Mg and Fe, on the other 
hand, could result in a sufficiently high T, but as we saw above they will not 
form sufficiently close to the star to drive a wind. 

At this point one might wonder about alternatives to the ballistic one-shock 
scenario for the trajectory prior to dust condensation. It might be possible that 
a matter element gets hit by several shocks at increasing distance from the star 
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if the fall-back is slow enough ('parachute regime', see above), gradually driving 
the gas out to a considerable distance (see, e.g., Bowen 1988). One could even 
think about replacing the first step(s) on the way with other mechanisms which 
on their own may not be sufficient to drive an outflow (just like the shock 
waves alone will hardly do the trick, at least not within observed constraints as 
known today). Several such possibilities have been discussed in the literature 
earlier, e.g., Alfven waves (e.g., Hartmann & MacGregor 1980), sound waves 
(e.g., Pijpers & Hearn 1989), or high thermal gas pressure due to heating by 
shock waves ('calorisphere', Willson & Bowen 1998). However, each of these 
alternatives faces the same question: will densities and dynamical timescales 
allow for efficient grain condensation at a certain distance from the star? 

2.3. Condensation efficiency and time scales 

Discussions about dust in atmospheres and winds of AGB stars often ignore 
the fact that a low enough temperature is only a necessary condition for the 
formation or survival of dust grains, and not a sufficient one. In dynamical at- 
mospheres and winds where timescales are set by pulsation and wind dynamics, 
the efficiency of dust formation is strongly dependent on prevailing gas densi- 
ties and abundances. As matter moves away from the star, the temperature 
decreases, which - at first - favors grain formation. At the same time, however, 
the growth of grains turns into a race against falling density which slows down 
the process, and more often than not leads to incomplete condensation. 

Simple estimates based on gas kinetics as presented in Gustafsson & Hofner 
(2004) demonstrate that the timescales for the growth of carbon grains close to 
the condensation radius (i.e. at about 2-3 stellar radii) are on the order of a year, 
i.e. comparable with the pulsation period of the star, and increasing outwards 
with falling density. This is in good agreement with detailed numerical models 
which tend to show mean degrees of condensation well below unity (typically 
0.3-0.5), and a rather limited zone of grain growth. 

Iron-poor silicates, which should have a condensation radius similar to 
amorphous carbon grains (see above), face the problem that the abundance 
of the key element Si is about an order of magnitude lower than that of C, in- 
creasing the estimate for the grain growth time by a corresponding factor (the 
timescale is inversely proportional to this abundance). This might result in a 
rather low dust formation efficiency, which, in combination with the low opacity 
of such grains, makes them a by-product of the wind, not a driver. 

Iron-rich silicates, on the other hand, which could in principle contribute 
to the total opacity, and consequently to driving the wind, will most likely not 
form in significant amounts, being handicapped by both the low abundance of 
Si and the much lower densities prevailing at distances corresponding to their 
large condensation radius. At such distances, the densities will be at least an 
order of magnitude lower than in the region where carbon or iron-poor silicates 
may form, adding another factor of ten to the grain growth timescale. 

3. Detailed models: status and trends 

The development of models for a particular phenomenon - in this case mass 
loss through dusty winds - often occurs in several steps: first ideas about basic 
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processes lead to order of magnitude estimates, followed by simple (analytical 
or numerical) models. If the basic principles seem sound the next step is an 
iterative improvement of numerical models by comparison with observations, 
leading eventually to detailed, reasonably realistic models which, finally, can be 
applied to study certain astrophysical phenomena in a wider context (e.g. the 
role of dusty winds in stellar and galactic evolution, to pick a not-so-random 
example) . 

Detailed time-dependent models for winds of AGB stars currently come 
in three major groups, namely two types of spherically symmetric frequency- 
dependent models including non-equilibrium dust formation for C-rich and O- 
rich chemistry, respectively, and 2D/3D models concerned with the effects of 
giant convection cells and structure formation on atmospheres and mass loss. 
These three types of models have reached different stages in their development, 
as will be discussed below0 

Early pioneering models in this field, exploring the general effects of pul- 
sations and dust for cool stellar winds (e.g. Wood 1979, Bowen 1988) make 
hardly a distinction between C-rich and O-rich stars, with the possible excep- 
tion of choosing appropriate values for certain input parameters. The turning 
point came with the inclusion of a detailed description of dust formation (e.g., 
Fleischer et al. 1992, Hofner & Dorfi 1997, Winters et al. 2000, Jeong et al. 
2003) in contrast to a parameterized description of the dust opacity, and/or 
the introduction of frequency-dependent radiative transfer, accounting for the 
complex opacities of molecules and dust (e.g., Hofner 1999, Woitke 2006a). 

The latest generation of models for atmospheres and winds of C-rich AGB 
stars by Hofner et al. (2003), combining a frequency-dependent treatment of 
radiative transfer with time-dependent hydrodynamics and a detailed descrip- 
tion of dust formation, compares well with various types of observations, such 
as low-resolution NIR spectra (Gautschy-Loidl et al. 2004) or profiles of CO 
vibration-rotation lines (Nowotny et al. 2005ab and this volume). With these 
non-grey dynamical models it is possible for the first time to simultaneously re- 
produce the time-dependent behavior of fundamental, first and second overtone 
vibration-rotation lines of CO, features originating in the outflow, dust forma- 
tion region, and pulsating atmosphere, respectively, probing the dynamics from 
the photosphere out into the wind. Recently, Mattsson et al. (2007) have ap- 
plied these models to investigate the formation of detached shells in connection 
with a He-shell flash. Currently, a large grid of dust-driven wind models for 
C-rich AGB stars and an accompanying library of variable synthetic spectra are 
being computed (see Mattsson et al., this volume). 

The development of similar detailed models for winds of O-rich AGB stars 
has been lagging behind the C-rich case, not the least due to a more complex 
scenario for dust formation. Jeong et al. (2003) presented wind models for 
M-type stars, combining a detailed description of dust formation with time- 



3 Other kinds of models, such as stationary wind models, neglecting pulsation and shocks, or 
pure dynamical atmosphere models without mass loss (e.g. Bessel et al. 1996, Scholz & Wood 
2000, Tej et al. 2003ab), are not included in this discussion. Note also that the recent models 
by Ferrarotti & Gail (2006) studying dust formation for a wide range of stellar parameters and 
chemical compositions are not wind models in the strict sense of the word since they have mass 
loss rates as an input parameter, not as a result. 
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dependent dynamics and grey radiative transfer. While the stellar parameters of 
these models are somewhat on the extreme side (high luminosities, low effective 
temperatures), the resulting wind characteristics are reasonably realistic. In 
view of the discussion in the previous section, it may be surprising that it is 
possible to drive winds with silicate grains. This is a direct consequence of the 
grey radiative transfer used in these models which corresponds to an effective 
value of p = (cf . Fig. [1]) , resulting in a rather small condensation radius for all 
types of grains, including iron-rich species. Frequency-dependent models for O- 
rich AGB stars by Woitke (2006a, and this volume) clearly demonstrate that the 
iron content of silicate grains has to be very low (large condensation radius for 
iron-rich silicates, combined with low condensation efficiency, see discussion in 
previous section) , and that the wind - even for quite extreme stellar parameters 
- will consequently not be driven by silicate grains, in contrast to previous 
expectations. 

During recent years, 2D/3D dynamical atmosphere and wind models for 
AGB stars have emerged, in addition to the spherically symmetric models dis- 
cussed above. The computational effort behind such models is considerable, and 
several simplifications have to be introduced in the description of physical pro- 
cesses. Nevertheless, investigating the effects of intrinsically three-dimensional 
phenomena like convection or flow instabilities on mass loss, seems a timely 
project in view of recent interferometric observations which indicate deviations 
from spherical symmetry. Woitke (2006b) presented 2D (axisymmetric) dust- 
driven wind models, including time-dependent dust formation and grey radiative 
transfer. He studied how instabilities in the dust formation process create intri- 
cate patterns in the circumstellar envelope, but without taking the pulsation of 
the central star into consideration. Freytag & Hofner (2003, 2007), on the other 
hand, investigate the effects of giant convection cells and of the resulting shock 
waves in the atmosphere on time-dependent dust formation in the framework of 
3D RHD 'star-in-a-box' models. The atmospheric patterns created by connec- 
tive motions are found to be reflected in the circumstellar dust distribution, due 
to the strong sensitivity of grain formation to temperatures and gas densities, 
as discussed above. 



4. Conclusions 

The well-known dichotomy between M-type and C-type AGB stars, as observed 
in molecular spectra, may have an even more drastic consequence for their mass 
loss mechanism. While advances in modelling, in particular the introduction 
of frequency-dependent radiative transfer in time-dependent dynamical models, 
have improved agreement between models and observations of C-rich AGB stars, 
the opposite seems to be true for the O-rich case. Recent models of M-type stars 
combining time-dependent dust formation with frequency-dependent radiative 
transfer demonstrate that silicate grains forming in such environments will be 
extremely iron-poor, resulting in too low opacities to drive a wind. 

These qualitative differences can be understood with simplified analytical 
considerations, as discussed here, causing serious doubt about the validity of the 
dust-driven wind scenario for M-type AGB stars, at least in its most simple form. 
The mystery is deepened by the fact that the observed wind characteristics for 
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both types of stars are rather similar (see, e.g., Olofsson 2004, Ramstedt et al. 
2006), which hints at a common mass loss mechanism. In addition, alternative 
scenarios discussed in the literature may have serious difficulties explaining the 
formation of considerable amounts of dust as a by-product in an outflow driven 
by a different force. Non-equilibrium dust condensation is very sensitive to the 
prevailing thermodynamical conditions, and restricted to a relatively narrow 
zone close to the star, putting strong constraints on potential driving forces. 

In this situation the role of observers should not be underestimated. Any 
observations which can narrow down the possible range of conditions in the wind 
acceleration zone are of great importance for solving this problem. 
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